


NON-INTERSECTING BROWNIAN MOTIONS LEAVING FROM AND 

GOING TO SEVERAL POINTS 

MARK ABLER, PIERRE VAN MOERBEKE, AND DIDIER VANDERSTICHELEN 

T-H Abstract. Consider n non-intersecting Brownian motions on R, depending on time t G [0, 1] , 

^^ with rrii particles forced to leave from at at time t — 0, 1 < i < q, and rij particles forced 

Cn to end up at bj at time t = 1, 1 < j < p. For arbitrary p and q, it is not known if 

^ the distribution of the positions of the non-intersecting Brownian particles at a given time 

Cl^ < t < 1, is the same as the joint distribution of the eigenvalues of a matrix ensemble. 

■^C This paper proves the existence, for general p and q, of a partial differential equation (PDE) 

satisfied by the log of the probability to find all the particles in a disjoint union of intervals 

E = U[^i[c2i-i, C2i] C R at a given time < t < 1. The variables are the coordinates of 

the starting and ending points of the particles, and the boundary points of the set E. The 

I I proof of the existence of such a PDE, using Virasoro constraints and the multicomponent 

Q^ KP hierarchy, is based on the method of elimination of the unwanted partials; that this is 

Q, possible is a miracle! Unfortunately we were unable to find its explicit expression. The case 

• p — q = 2 will be discussed in the last section. 
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2 m. adler, p. van moerbeke, and d. vanderstichelen 

1. Introduction 

Consider N independent Brownian motions on the real line, starting at time t = and 
ending at time t = 1 at prescribed positions, and conditioned not to intersect during the 
time interval ]0, 1[. In the particular case when all the Brownian motions start at the same 
position at t = and end at the same position at t = 1, let's say the origin, the positions 
of the Brownian particles at an intermediate time < t < 1 have the same distribution, 
after a space-time transformation, as the eigenvalues of a randomly chosen matrix of the 
GUE ensemble. It is also the distribution of A^ Dyson Brownian motions on the real line. 
This process, discovered by Dyson [15j, describes the motion in time of the eigenvalues of a 
N X N Hermitian matrix whose real and imaginary parts of the entries perform independent 
Ornstein-Uhlenbeck-processes, with an initial distribution given by the invariant measure 
of the process. In the case of one starting position and two (or more) ending positions, 
the positions of the Brownian particles at an intermediate time < t < 1 have the same 
distribution, after a space-time transformation, as the eigenvalues of a randomly chosen matrix 
of the Gaussian Hermitian ensemble with external source. The relationship between non- 
intersecting Brownian motions and matrix models has been developped by Johansson |16| . 
and for the Gaussian Hermitian matrix ensemble with external source by Aptekarev-Bleher- 
Kuijlaars |9j. See also Adler-Delepine-van Moerbeke [Ij and Katori-Tanemura |18^ I19j for a 
detailed description of the relationship between Dyson Brownian motions, non-intersecting 
Brownian motions and Gaussian Hermitian matrix ensembles. In particular, in [18] both 
stochastic processes are obtained as scaling limits of the vicious walkers model. In the two 
particular cases cited (i.e. non-intersecting Brownian motions with one starting position and 
one or several ending positions), the relationship between non-intersecting Brownian motions 
and Hermitian matrix models has led to a deeper comprehension of the diffusion problems. In 
both cases, partial differential equations (PDE) for the finite N diffusions have been obtained 
(see [HEIIH]). For large N, upon taking appropriate scaling limits, different processes appear 
describing the transition probabilities of critical infinite dimensional diffusions, like the Airy 
process, the Sine process and the Dyson process (see |71[T9l[2j) for one starting and ending 
position, and for two or more ending positions the Pearcey process (see O [22]), the Airy 
process with k outliers (see [1]), etc. The approach of Adler-van Moerbeke is to take scaling 
limits of the PDE's describing the finite N processes to obtain PDE's for the transition 
probabilities of the critical infinite dimensional diffusions. 

Consider now N non-intersecting Brownian motions on the real line starting at time t = 
at q and ending at time t = 1 at p prescribed positions, with p,q > 2. It is not known if 
there exists a matrix ensemble whose joint eigenvalue distribution describes the distribution 
of the positions of the Brownian motions at an intermediate time < t < 1. For p = q = 2 
this problem has first been studied by Daems-Kuijlaars |10j and Daems-Kuijlaars-Veys |llj . 
In these papers, the authors consider A''/2 particles going from a to 6, and A^/2 particles 
going from —a to —b. They show that the correlation functions of the positions of the 
non-intersecting Brownian motions have a determinantal form, with a kernel that can be 
expressed in terms of mixed multiple Hermite polynomials. They analyze the kernel in the 
large A^ limit, for a small separation of the starting and ending positions (i.e. when the 
product ab is sufficiently small) , and find the limiting mean density of particles is supported 
by one or two intervals. Taking usual scaling limits of the kernel in the bulk and near the 
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edges they find the Sine and the Airy kernel. For large separation of the starting and ending 
positions, those results have been extended by Delvaux-Kuijlaars fT2]. In [2], Adler- Ferrari- 
van Moerbeke study a similar situation, but with an asymmetric number of paths in the left 
and right starting and ending positions. Recently, Adler-Ferrari-van Moerbeke [3j and also 
Delvaux-Kuijlaars-Zhang [T3j (see also [13]) analyzed the large A^-limit in a critical regime 
where the paths fill two tangent ellipses in the time-space plane. Using an appropriate double 
scaling limit, they prove the existence of a new process describing the diffusion of the particles 
near the point of tangency. 

It seems to be a highly non-trivial problem to obtain concrete results about the processes 
describing the critical infinite dimensional diffusions, obtained as limiting situations of the 
problem of N non intersecting Brownian motions on the real line starting at q and ending at p 
prescribed positions, with p,q > 2. The aim of this paper is to provide a better understanding 
of the finite A^ diffusion for two or more starting and ending positions. We consider A^ non- 
intersecting Brownian motions xi{t), . . . ,xi\i{t) on M, starting at time t = in q different 
points, and arriving in i = 1 in p different points. If P^^'"'jj''(all Xi{t) £ E) denotes the 
probability to find all the particles in a set E at an intermediate time < t < 1, we prove the 
following theorem. 

Theorem 1.1. For each value of the parameters p > 1 and q > 1, let K* he the smallest 
positive integer such that 

(x^ - 3x + 4)(K*)2 + (-x^ + 3x + A)K* - 2x{x^ - 2x - 1) > 0, 

with X = p + q. Let E he a finite union of intervals. Under the assumptions ai + • • • + a^ = 
and hi + ■ ■ ■ + hp = Q, the function logFj^^''"'^' (a// Xi{t) G E^ satisfies a nonlinear PDE of 
order K* + 3 or less, the variahles heing the coordinates of the endpoints of the set E, and the 
coordinates of ai, . . . ,aq and hi, ... ,bp. 

For example, for 4 < x < 8, the value of K* in this theorem is given in the following table : 



X 


4 


5 


6 


7 


8 


K* 


3 


4 


5 


5 


5 



The proof of Theorem 1 1 . 1 1 will be given in section 5, and is based on the use of a particular 
integrable hierarchy, and Virasoro constraints. The use of these methods is suggested by the 
fact that the probability P^^''"'^' (all Xi{t) G E^ has different descriptions (see section 2): 

(1) It can be written, after making a space and time transformation, as a block moment 
matrix 



tZZ{^^^-^it)^E)=^ 



det 



x™^,(x) y>,(y) 



0<m<m,i — 1 



l<i<q 
l<j<P 



(1.1) 



where ipii^) = e""*^, 'fj{y) = e ^^, and the following inner product 



x™V^(x) y>,(y) 



^m+n^{ai+bj)x^- 



2 dx. 



The ~'s indicate that a space-time transformation has been performed (see section 2, 



formula (2.2)). 
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(2) It can be written as a sum of multiple integrals 

a^bf^ 4=1 «=i 

nig 

X ... X IAmg(3;miH |-mg_i+l5 • • • ) S^miH f-rrtq j J_ J_ ''/-'ij (.^Jmi 



Za 



H |-m,-i+J, 



i=l 



"■1 



X • • • X 



i=l 

(^^np(a^o-(ni+---+r!,p_i+l)) • • • ) 2;o-(ni+---+np)) J_ J_ '/'p(2^cr(ni+--+np_i+j)) J ) (1-2) 



i=l 



where A„ is the Vandermonde determinant, and S^ is the group of permutations of 
A^ elements. 

As shown in Adler-van Moerbeke-Vanhaecke |8] , the determinants of block moment matrices 
deformed in an appropriate way satisfy integrable hierarchies. Concretely, the determinant 



(1.1) is deformed by adding exponentials containing additional families of time variables. 



one family for each weight function ^pi and ■i/'j) or equivalently, (1.2) is deformed by adding 



exponentials containing additional families of time variables, one family for each Vandermonde 
determinant. The determinants of the deformed block- moment matrices (1.1) are then tau 
functions for the multi-component KP hierarchy. The multi-component KP hierarchy is a 
very general hierarchy of integrable equations, describing the time-evolution of matrix- valued 
pseudo-differential operators, depending on several families of time variables. These operators 
can be expressed in terms of so-called tau- functions, which encode the whole hierarchy. As a 
consequence, the determinants of the deformed block moment matrices satisfy some nonlinear 
PDE's. This is developped in section 3. 

When p = 1 or q = 1, it is easy to see that all the terms in (1.2) are equal to each other. 



and the sum is simply A^! times a A^-tuple integral over E. This unique integral corresponds 
to the joint eigenvalue probability of the Gaussian Hermitian ensemble with external source. 
It is a well known fact that matrix integrals deformed in an appropriate way satisfy Virasoro 
constraints (see [1]). These constraints are linear PDE's satisfied by the deformed integrals. 



involving a time part and a boundary part. Although we do not know if (1.2) for general p 



and q corresponds to (the reduction to polar coordinates of) a matrix integral, we show that 



each term in (1.2) separately satisfies Virasoro constraints. As a surprise, it appears that all 



the terms satisfy the same Virasoro constraints, and hence, by linearity, it follows that (1.2) 
satisfies Virasoro constraints. This is developped in section 4. 

Following a method developped by Adler-Shiota-van Moerbeke (see for example [4,j and 
[5]), Virasoro constraints with time and boundary parts can be used to eliminate all the 
partial derivatives with respect to the added time variables in the non-linear PDE's from 
the integrable hierarchy, and hence to obtain a non-linear PDE with respect to the variables 
of the unperturbed problem. The complexity of the problem studied in this paper does not 
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enable one to perform concretely this elimination process and to obtain an explicit formula for 
arbitrary values p,q > 2. It is a priori not even obvious at all that it converges to a PDE after 



a finite number of steps! In Theorem 1.1 we prove, however, using a simple combinatorial 
argument, that it indeed does, and this for general p and q. We would like to emphasize that 
the existence of a PDE satisfied by Pj^'"'(^''(all Xi{t) G E^ is not obvious at all. Our proof 
rests on two surprising facts, the first being that the perturbed problem satisfies Virasoro 
constraints, and the second that the elimination process converges after a finite number of 
steps. 

In the last section, we consider the case of non-intersecting Brownian motions with two 
starting and two ending positions. In this particular case, we improve considerably the general 
method given in the proof of Theorem |1.1| to obtain a PDE. 



2. NON-INTERSECTING BrOWNIAN MOTIONS WITH q STARTING POINTS AND p ENDING 

POINTS 

Consider N non-intersecting Brownian motions xi{t), 2^2 (t), . . . , xj^{t) in M, leaving from 
distinct points ai < 02 < • • • < Qn and forced to end up at distinct points /3i < /32 < • • • < /Sat . 
Prom the Karlin-McGregor formula |17j . we know that the probability that all Xi{t) belong 
to a set ii^ C M at a given time < t < 1 can be expressed in terms of the Gaussian transition 
probability 

, , 1 i^-y)^ 

p{t,x,y) = —=e t , 
vvrt 



in the following way 

P^(an Xi{t) G E) := P( all Xi{t) G E 



(2:1(0), . . . , xn{^)) = (ai, • • • , an) 
(xi(l),...,x^(l)) = (/3i,...,/3^) 



Zn 



N 



1 

Zn Je'^ 



/ det [p{t, a,, Xj)] ^^. det [p{l - t, Xi, /3j)] ^< . TT dx 



where Z^ and Z^ are normalizing factors. In particular, if 

(ai,...,a7v) = (ai,ai,. .. ,01,02,02,... ,02,. .. ,ag,aq,... ,ag), oi < 02 < • • • < Q 

mi m2 m, 

(/3i,...,/3Ar) = {bi,bi,...,bi,b2,b2,...,b2,...,bp,bp,...,bp), 61 < 62 < • • • < 6p, 



"1 
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with XlLi "-i = Ya=i bi = and Ylj^i rrii = Yl^^i "« = ^) then we have 

<".:;; (all x.(i)Gi5;) 



/ 



all Xi{t) G E 



V 



(xi(0), . . . , xn{0)) = {ai,...,ai,. . .,aq,. ..,aq 



)\ 



nil 



(a:i(l),...,XAr(l)) = {bi,...,bi,...,bp,...,bp) 
^ ^ ^ ^ ^ 



ni 



lim 

Ql,...,Qmi — ^ai 



$(all x,{t) £ E) 



■/3 



/3l,...,/3„-^^6l 



Consequently, we obtain 



1 



"Zzz (all x,{t)eE)=^ I n ^"^ '^^^ 



Af 2 



/ 



X det 



^j^ * )o<i<mi-l 



l<j<N I 



. det 






2bpXj 

^j^ ^ * )0<i<np~l I 



We have, using the change of variables Xj 



%^y^, i<^<iv, 



t:.:,t(aii^.we^) = ^p,. 



'^-V^«Vi^' 



(2.1) 



with the normalized problem being 



Pp,g(^;a,6):= 



Z, 



P,? JE'^ 



N 2 

J_J_e^^ dxjj det [^i(a;j)] i<._^.<^ det [(pi{xj)] 



1<«J<A^' 



j=l 



(2.2) 



where Zp^q is a normalizing factor, and where we have introduced the following notation 



(Vi(x),...,V'^(x)) := e'^l^xe 



,a.ix ^^aix 



.. . ,X 



rrii-l aix a2X ^^a^x 



xe 



j.m2-l^a2X 



f^^qX rpfp^qX 



rf-mq-l^agX 



{ifl{x),...,ifN(.x)) 



pbix ^J>\x ^rii-l bix b2X ^Mx 



n2-lpfe2a; „bpX ^p^P^ ^"p-lpfep^ 
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In the following proposition, we consider a general situation, of which (2.2) is a special case 
by setting V{x) = ^ , ipi{x) = e"*^, I < i < q, and (pi{x) = e^*^, 1 < i < p. 

Proposition 2.1. Given an arbitrary potential V{x) and arbitrary functions ?/'i(^)) • • • ^i^qi^) 
and (pi{x), . . . , (Pp{x), define (N = mi + • • • + nriq = ui + . . . Up) 



(Vii(x),...,-0Ar(x)) ■.=[il)i{x),x%lji{x),...,x'^^ Vi(a;),V'2(x),xV'2(a;), 

• • • , X™2- V2(X), . . . , ^q(x), XV,(X), . . . , X'"'- V,(^)) , 
((^l(x),. . .,(Pn{x)) ■.= Upi{x),Xipi{x),. . .,X^^^'^ipi{x),ip2{x),XLp2{x),... , X^'^~^ ip2{3 

tpp{x),xipp{x), ..., x^-p'^ifpixU . 



We have 



I' ^ 



(A^!)2det 



Tpi{x)if>j{x)e ^'^^^dx 



J l<i,j<N 



N \ / N 

mi,m2,...,mgj \ni,n2,. 

rriq 



X . . . X 



/ (T\e~^(^^Ux^(ArnAx^'^)T\4^ 

1=1 

X] (-1)" [(^«i(^<x(l)' • • • ' 2;a(ni)) n V'l(a^<7(i)) 

i=\ 

(^^np(a^cr(ni+---+np_i+l)) • • • ) 2;o-(ni-|— -|-np)) J_ J_ '/'p(2;cr(ni+---+np_i+i)) 



X • • • X ( A„„ ix 



i=X 



(2.3) 

where x^^'^ = (xi, X2, . . . , x^J, •■ ., x^^) = (xmi+-+m,_i+i, • • • , Xmi+'-fm,); andA„(xi,...,x„) 
det[x*^ ]i<i,j<n lis ^/ie Vandermonde determinant. 



Proof. Let P{E;p,q) be left hand side in (2.3) 



P{E;p,q) ■=N\J^ ([Je-^^^') dx^) det [Vii(a:,)],<.^.<^ det [^iixj] 



'i<i,j<N' 



The first identity in (2.3) is a consequence of applying the following standard identity 
and distributing the integration over the different columns. 
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We prove now the second identity in (2.3). Working out the determinant det [^i(xj)] ^^^ . 
we obtain 



<N 



P{E;p,q) =N\Y, (-^T / ^ (H e"""^""^ ^^i) {Il^^iM^))) ^et ['P^{x,)]^^ 



In each term of this summation, we make the change of variables Xi = Ha'^u), 1 < i < N. We 
have 



» iV N 

P{E;p,q) = {N\f / [He-'^'^y^Um) [l[My^)) det [^,{y,] 



since det [<^i{yj)]i^ij<p^ = (-l)'^det ['Piixj)]-^^^.^.^^. Now take cti G S^i, cr2 S Sm2, ••■, 
aq £ Smq arbitrarily and define the permutation a := cri x c72 x • • • x cjg G Sn- Consider the 
following change of variables y ^- z defined by a: 

(yi,- ■■,ymi) = iZai{l),-- ■ ,^ai(mi))) 
(,2/mi+l) • • • 7 ymi+m2 j = \Zmi+a-2{l)'> • • ' ' '^mi+o-2(m2))' 



[ymi-\ hm,_i+l5 • • • 5 ymiH hm, j — [Zmi-\ hmq_i+(Tq(l) i • • • ) ^miH [-mq_i+o-q(m,) j- 

This change of variables leaves the integral unchanged. Consequently, if we sum over all the 
permutations ui x (T2 x • • • x cjg G Smi x • • • x 5^, and divide by mi\m2l ■ ■ ■ ruql, we have by 
the definition of ipi and An{z) 

mi, 7712, ... , nigj Jen \fj^ J\ f-^ 

1712 "^9 

x(^A„^.,{z^'^^)'[lMzm^+i)) X ••• X (^A™,(z('?))JJV'g(^mi+...+m,„i+i))det[(^i(zj)]i<._^.<^. 

(2.4) 
We develop the determinant det [<^i(-2j)] ]^<j -^^ in this expression 



N 



det [<^i(2j)]i<i,,<^ = X] (-i)"n'^*(Mo) 



crSS'jv *=1 



ni n2 

rip 

4=1 
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Fix o"! G Sni, • • • , <5"p G Snp and, for a given permutation a G Sj\f, let it G Sat be such that 
a = a o (ai x cj2 x • • • x ap), but when 5" runs over Sn, then so does a. Consequently we have 



det [(fiAz.i 



Z^r.;^,(i\]zl\ ,,. IX.., 






"1 






j=l 



We substitute this expression in equation (2.4). For each a G Sj\f we further sum over 



o"! G 5„^, o"2 G 5^2, • • • , fTp G S'np and so must divide by ni\n2l ■ ■ ■ Upl as we have overcounted. 
We then obtain the second equality in (2.3). This ends the proof. D 



3. An INTEGRABLE deformation and {p + g)-COMPONENT KP 

The connection of the problem of non-intersecting brownian motions on M with the multi- 
component KP hierarchy is explained in [8| . The main ideas are being sketched in this section. 

3.1. The (p + g)-component KP hierarchy. Define two sets of weights 

V'i(x),...,'0g(x), and ipi{y),...,ipp{y), with x,yGM, 



(4' 



(") Ac) 



) •'2 ' 



. . . ), 1 < a < g, and 



^-W..k 



and deformed weights depending on time parameters s^"^ 
tW = {tf\tf\ . . . ) , 1 < /3 < p, denoted by 

iIj~'{x) := 'il;a{x)e-^'^-^'k" ^\ and c^^(y) := 99/3(^)6^"=!**'''^'. 
For each set of integers 

rh = {mi, . . . ,mq), n = {ni, . . . ,np), with |m| = |n|, 

\m\ = Yl'i=i''^i^ 1^1 ~ X]?=i%' consider the determinant of the moment matrix T^^ of size 
\m\ = \n\, composed of pg blocks of sizes miUj. The moments are taken with regard to a (not 
necessarily symmetric) inner product (• | •) 






detT^ 
/ 



det 



{x^^-^'{x)\y^^{{y)) 



0<i<m\ 
0<j<ni 



0<j<np 



\ 



0<j<711 



( (a^^V'/W I yVp(y)> )o<i<m, 

0<j<np ) 



(3.1) 
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Theorem 3.1 (Adler, van Moerbeke, Vanhaecke |8]). The determinants of the block matrices 
Tffifi satisfy the bilinear relation^Q 

for all m, n, fh* , n* such that \fh*\ = |n*| + 1 and \rh\ = \n\ — 1, and all s, t, s* ,t* G C°° , and 
where 

a p 

aa{rh) = ^ {ma' - m^,), and (Tp{n) = ^ {upi -n*^,). 

a'=l /3'=1 

These identities define the {p + (7)-component KP hierarchy, as described by Ueno and 
Takasaki EOl. 



Define the Hirota symbol between functions / = f{ti,t2, ■ ■ ■) and g = g{ti,t2, ■ ■ ■), given a 
polynomial p{ti,t2, ■ ■ ■), namely 



y=o 



Kdt;^dt2--)f°''-=Kdy-i'dy-2--)^^'^'^'^'~'^ 
This operation extends readily to the case where p{ti ,t2, ■ ■ ■) isa Taylor series in ii , t2 , • • • • We 
also need the elementary Schur polynomials si, which are defined by e^fc=i ^^^ := X^fcLo Sfc(0-^'^) 
for / > 0, and si{t) = for / < 0. Moreover, set 

(91915 
^dh' 2dt^' Sdts' ^ 

With these notations, computing the residues about z = oo in the contour integrals above, 
the functions r^^, with \m\ = \n\, are found to satisfy the following PDE's : 



fa\ f d I d 1 d \ 



tI 



52 



rhn m\ tpi\ ^^S ''mn — Sj-\-25ppi\Pt(ll)Prn,n+ep-eiJ, ° ''"m.n-e/j+e^/ ; 



2 



'ds^^-^ds\^ 



,, lOgTf^fi — Sj^26^^i\'^s(°')rrh-ea+e^i,n ° 'Tfh+ea-e^i ,ni 



r^^ — 77^ logr^ri = -Sj{d^iP))Trfi+g^^fi+ef, ° T« 



dtl'dsWi 



Introduce the notation [z\ = {z, ^, ^, . . .) for 2 G C, and let ei — (1, 0, 0, . . . ), 62 = (0, 1, 0, ... ), .... 
Only shifted times will be made explicit in the functions r^^. The integrals are contour integrals along a small 
circle about 00, with formal Laurent series as integrand. 
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3.2. An integrable deformation of the joint probability density function for the 
problem of non-intersecting Bro^vnian motions. We will now deform Pp^q{E;a,b) de- 



fined in (2.2) by adding extra time variables 



.(1) _ (,W fW \ A2) _ (A2) .(2) N Jp) _ (Ap) Ap) ^ 

b — yb^ , ^2 1 ■ ■ ■ J 1 ^ — \^l ) ''2 T ■ ■ ■ j 1 ■ ■ ■ ) '' ~ v''l ' ''2 t ■ ■ ■ ) 1 
Jl) _ fJl) Jl) ^ J2) _ / (2) (2) N Iq) _ , iq) iq) ^ 

and auxiliary variables 

(ai,...,ag), (/?!,.. .,/3p), 
such that ^1^1 Oii = Yl^j=i l^j = 0- First set 

yjj^ [xj := e -^j^i j ^ 1 < ^ < g, 

^2 I Y-oo ,(i)^ 






We define 



T"? ^(t, s] a, 13; a, b) 
Pp,g{E;a,b;{t,s),{a,p)) = -!^^ -, (3.3) 



with 



rifi{t,s;a,l3]a,b) ■= —I (^JJe^dxij det [^^['{xj)]^^^^^]^ det [(^-(xj) 



i=l 

(3.4) 



where {t,s) = (t«, . . . , t(P); s«, . . . , s^^)), (a,/3) = (ai, . . . , a,_i; /3i, . . . ,/3p_i), (a,6) 
(ai, . . . , aq_i; 6i, . . . , 6p_i), and 

(V;r(x),...,Vi^^(x)):=(v^r(x),x^r(x),...,x"^-Vr(x),..., 

V'g""(x),xV^-^(x), . . . ,x™«-V<^'(a;)), 
(c^*i(x), . . . , c^*^(x)) :=(^*i(a;), Vi(a;), • • • , x^^^VK^;), • • • , 

99^(x),x(^^(x),...,x"''-Vp(a;)). 
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Observe that Pp^q{E;a,b) = Pp^q[E;a,b;{t,s),{a,/3))\^, where L = {{t,s) 
0}. By virtue of proposition |2. 1| we have 



(0,0), a = /? 



'^m.nl*' s; «' /5; «' &) = det 



ip- "" {x)(p^j{x)e 2 dx 



i<i,j<N 



det 



0<i<»ni 
0<j<ni 



/gX*+^VqV*ie "2 c?x)o< 



l<mq 

0<j<ni 



0<J<71p 



0<i<np / 



(3.5) 



where for simphcity we have left out the dependence of ip^ ^ and ijj^- on x. We have also 



TiJt,s;a,P;a,b) 



miH hruq 



mi 



An,i{x^^'^)Ylipi'{xi)e 2' dx, 



X . . . X 



i=l 



(a^,(x('^)) J] vr(^i)e^^^' 



i=miA hmq_i + l 



m 



X • • • X 



(t£5jv 4=1 

niH hrip 

^np(^o-(niH hrip-i+l); • ■ ■ ' ^cr(ni+---+np)) J_ J_ ¥'p(3^o-(j)) 

i=niH hnp-i + 1 



(3.6) 



fj^f{z)g{z)e ^ /^^z, with 



The determinant of the moment matrix (3.1) with regard to the inner product {f,g) 



ipi{x) := e' 



aiX+aiX" 



l<i<q, 



J>jX+l3jX- 



^,{x):=e'^-^^^-' , \<2<V. 



is the same as the determinant O.Sh. Therefore, by virtue of Theorem 3.1 , t^ ^{t, s; a, (3; a, b) 
satisfies the {p + g)-component KP hierarchy. 
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Corollary 3.2. The function T^f^{t,s;a, /3;a,b) satisfies the following identities, I < k,k' < 
q,l< l,r <p, k^k', l^V, 

jp 9 1 E rp d 1 E 

ft T^ ^ , ^ ^ fiiWaA'') rn,n f) T^ ^ , ^ ^ a {ft)Q (fc') '''■'■' rh.n 



Oil m,n-ei+ei, Qt^Og^C) '^'■'^ ' rn,n OS-^ rn+ek-e^,,n g^^'^'as''''^ '-'■'■' rn,n 

C 9 In T^ T7 9 ]„ _-E 

f) r-5',-. -.^ fl/Wflct*) m,n f) T-5' - -, - ;5^(*);5+(0 "^ 'm,n 

aAl) T^ d^ inr^ ' f}M t? ^ ^ ^ 9^ ]„ ^^ ' ^ ^ 

Proof. We shall only give the proof of the first identity. The two first elementary Schur 
polynomials are given by 

So{xi,X2,...) = 1, si(xi,X2,...) =a;i. 



Consequently, the first equation in (3.2) with j = and / 7^ /' gives 



:)2 



a2 
/E\2t/,£;_,^N£; e _ E E ("i ii\ 

Vmn) n\ (;/■) ■^*-'§ ''"mn ^^K^tO-) )'^m,,n+ei—eii ° ''"m,n— ej+e";/ ^m/n+ej— e";/ ^m,n— ej+e;/ ' V"-"J 

while for j = 1 and I ^ I' \t gives 
;^2 

\'mn) fn (;/■) 6 ' mn ''1 v'^tW J 'm,n+e;— e,/ ''m,n—ei+e,i 

at^'oti ' 

= T^ ^ E _ E ^ E /o q^ 



Taking the ratio of (3.9) and (3.8) yields the first formula of Corollary 3.2, The other identities 



are obtained in a similar way. D 

4. ViRASORO CONSTRAINTS 
Let us introduce the following differential operators 

Otm 

i+j=m -' i>l i+j=—m 

/, m + l\/ (9 , , \ k(k + l) ^ 
+ (^ + ^) (^ + (-m)t_„) + A_^ ,^,. 

Those operators satisfy the Heisenberg and Virasoro algebra respectively 
and interact as follows 

We have the following lemma, proven by Adler and van Moerbeke [6]. 
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Lemma 4.1 (Adler-van Moerbeke [6]). Given p{z) = e~ '^•', with 

p'{z) _ ^, _ g{z) _ EZo^iz' 



P{z) 



f{z) T^ZoP'iZ' 



the integrand 



N 

dlN{z;t) := An{z) JJ (e^»=i*'^^/9(zfc)dzfc) 



fc=i 



satisfies the variational formula 
d 



de 



dlN{zi ^z^ + efizi)z'l+';t) _=Y1 (w4+/,JvW " ^aii+i,jvWV^^iv(^;t), 



e=0 



/=0 



^(1) 



for each k > —1. The contribution of the factor Yl^—i dzi in this equation is 

oo 

Y,Piil + k + l)/^li,NdlN{z;t). 



1=0 



We define, for a given permutation a & Sn, tfie integrands 



mi 



dlZ^f,{x;{t,s)) =(A„,(x«)nV'i \x,)e^dx 



X . . . X 



j=l 



miH \-mq 2 

(A™,(x('?)) n ^p-'"'\x,)e^ dxi) 

i=miA hm,_i+l 

rii 

X (^ni (a;<^(i), • • • , 2;^{ni)) H '^i^" (^^tw)) x • • • x 

i=l 

niH hrip 

(A„p(2;^(ni+...+np_i+l),---,a;„(ni+...+np)) JJ V'p''' (a^^w)) • (^-l) 

i=niH hrip-i+l 

We are looking for a variational equation for 

dl^,n{xi ^Xi + ex^+^- {t, s)) . 
We have the following lemma. 



Lemma 4.2. The integrand dIZ ^(^x; {t, s)^ as defined in {4-1), satisfies the following varia- 
tional equation for each a G Sn and k > —1 






^ m ni'^i '~~^ "^i ' ^"^i i V ' ^ 



e=0 
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with 



V 






(4.2) 



j=l 



Proof. By the Leibniz rule, applying Lemma 4.1 to each factor in (4.1) and adding all these 
contributions yields (4.2). D 



As the variational formula in Lemma 4.1 is independent of the labeling of the variabl es in 
the integrand, it is a trivial but very important fact that the operator V™'" as defined in (4.2) 
is independent of the choice of o" G 5"^. As a consequence, we have the following theorem. 



Theorem 4.3. The function T^^{t,s) as defined in (3.4) satisfies the following Virasoro 
constraints 



^fe ^rn,n 



with 



ym,n E 
" k m,n^ 



2r 



k> -1, 



(4.3) 






dci 



forE = [Jl^-^[c2i-uC2i] CM. 



Proof. By virtue of formula (3.6) expressing T^f.{t,s;a,f3;a,b) as a A^-uple integral over E, 
we have 



^m,n(^''S;">/3;a,&) 



UUm^lllU^^ 



Y^{-ir4'Zit,s;a,P;a,b), 



cre5iv 



where T-.'Z{t, s; a, /3; a, b) is defined by 



E,a 



Trn'rii^^s;a,^;a,b) 



di^Ax;{t,s)), 



with dIZ_^(x; {t,s)) as in (4.1). For a fixed permutation a E Sn and k > —1, we apply the 
change of variables Xj i— )• Xj + ex^ , 1 < "i < -/V, given in lemma 4.2, in the integral defining 
T^'Z(t, s; a, P; a, b). This change of variables leaves the integral invariant, but induces a change 
of limits of integration, given by the inverse map 

Ci^Ci-ec,^+^ + 0(e2), 1 < i < 2r, 

for e small enough. Consequently, diff^erentiating the result with respect to e and evaluating 
it at e = 0, using the fundamental theorem of integral calculus together with Lemma 4.2, we 
obtain 



''■ m,n k, m,n^ 



k> -1, 



(4.4) 
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with 



St 



2r 

i=l 



,fc+l 



d_ 
dci 



As noticed earlier, the operator V™'" does not depend on cr G Sn- Consequently, summing 



(4.4) over a G Sn and dividing by ni'=i "^i- Tlfci "-i'' '^^ obtain 

2 S _^,m,n,£ 



V 



/c> -1. 



This concludes the proof. 



D 



When specializing the differential equations (4.3) to k 
tau- function t^ . 



-1 and /c = 0, we find that the 



S_ir 



satisfies respectively 
d 



Bor 




+ 2^A 



dr 



dt 



(0 
1 



(4.5) 



log 



The Virasoro constraints (|4.3|) play a very crucial role in finding a PDE for the function 



bi, 



[all Xi{t) G E) in the variables ai, . . . ,ag, bi, . . . ,bp and the endpoints of the set E. 



In the next section, we will prove the existence of a PDE for the logarithm of the normalized 
problem Pp^q{E;a,b) defined in (2.2), and deduce Theorem |1.1| from it. The normalized 
problem is related to the function t^ ^{t, s; a, (3; a, b) on the locus L = {{t,s) = 0, (a,/3) = 
0} through formula (3.3). As we have seen, this function is a tau- function of the {p + q)- 
component KP hierarchy and thus satisfies the PDE's (3.9). As the Virasoro constraints 



involve derivatives with respect to the endpoints of the set E, as well as derivatives with 
respect to the time variables (t, s), we will prove that they can be used to eliminate all the 



derivatives with respect to the time variables in (3.9) on the locus L. The proof proceeds in 



two main steps. In the first step, the Virasoro constraints, together with the linear conditions 
imposed on Oj, Oj, 5j, /3j, 1 < i < q and I < j < p 

q p q 



=1 



=1 



i=l 



E/3^ 



(4.6) 



are used to express on the locus % = {{t,s) = 0} all the derivatives with respect to the time 



variables in (3.9) in terms of derivatives with respect to the auxiliary variables ai, . . . ,aq and 
/3i,...,/3p. In the second step, using a combinatorial argument, it will be shown that, on 
the locus L, all these derivatives with respect to ai, . . . , a^ and /?i, . . . , /3p can be eliminated. 
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Both steps will be performed in the next section. We end this section with some consequences 
of Theorem 14. 31 

From the linear conditions (4.6) it follows that the function T^^{t,s;a, (3;a,b) as defined 



in (3.4) satisfies the following equations 



1=1 CSj 1=1 CCj 



0, i>l, 



and 



dai 



dr dr 



Si) 



dr dr dr 



(<?)' 



dr _ dr dr 

dr _ dr dr 



, 1 < i<J3-l, 
, l<i<q-l, 

2 

, l<i<p-l. 



(4.7) 

(4.8) 

(4.9) 

(4.10) 

(4.11) 



From these equations, we deduce two families of identities. Let / := log rj^ ^(t, s; a, /?; a, 6). 

df 



Firstly, using equation (4.8) we have 

<? 

1=1 (y^i 

since X]i'=i ^i — 0) ^^"^ similarly 



q-l 

1=1 



dai'' 



(4.12) 



>f df 
1=1 ds\' 



g-l 



1=1 






p-i 



E 



ar 



1=1 ds^ 



df_ 



V-, df ^'df 

1=1 C^i 1=1 ' 

g-l 



df 



p-1 



EA£fe = EA 



1=1 (^^1 



1=1 



dl 

dbr 



1=1 



dai' 



S-o df x^^df 

1=1 0^2 1=1 ' 



(4.13) 



Secondly, using equation (4.8) we have 



9-1 



^JUUsr^df , df 

Z^ r, (i) 2^ 



:t asf ^ ■^ dai Qgi^ 



(<?)' 



and thus 



df ^i' df iV.df 



Using again equation (4.8), we obtain 






(4.14) 



(4.15) 
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Equations (4.14) and (4.15) can be summarized as follows 



df 1^9/ 1^5/ 



df 



ds[^^ 1 ^ dsf Q ^ 9ai 
Similarly, we have 



E^-^E^-c-Wi-. i<i< 



(4.16) 



df i^ df 1 ^4 a/ 



df 



Eii)--Ei:-(i-^.p)i:' 1^^-^^' 



df Is^ df 1^9/ 



df 



dsi'^ 1 U ds 



1=1 



5q/ 



9a,- 



^ _ 1^ J^_ li^ 9/ _ _ 5/ 



(4.17) 



Substituting relations (4.12),( |4.13[ ), ( |4.16[ ), ( |4.17[ ) in the Virasoro constraints ( |4.5[ ), we get 

/ 1 



•^■Sl ^ i>2 \ i=l 



IS- 



(0 ^ 



9s 



(0 






P 



1 / 

+ - ( ^ aimi + ^ 6;ni ) , I < j < Q-, 

1=1 
/ 9 



^./ 



"^ 1=1 
df , 1 






E E 



iS 



JO 



1 1=1 (^'^i-iJ ^ 1=1 1=1 



1/ '' ^ \ 



P 1=1 
df , 1 



1=1 
q 



^./ = -fe + ^E(E-S"i^+E"l"' 



K 



'2 



'-'^o ^ i>l \ 1=1 



Bjf 



df , 1 



dsY 



1=1 



dt 



ii^i^ + T '-'■-'' 



d 



d 



4^-) 



^^o <■ i>i \ 1=1 



EIE-l"^ + E«i"T^i/+ 



9s 



/=i 



dt) 



(0 



ir 

! 

P 



1 < i < p, 



(4.18) 



where 



g-l 



z=i 
P-i 



'"' a ^9 



^^--(^-^-)^nr-+^^-KE-^+EA,, 



'-' d ^9 



^^--(^-^-)|;%-h"^|-H^"'a^"^^^. 



1 < i < g, 
1 < J < p, 



/=i 



Z=l 



i=l 
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and 



F) ^( ^'^'^ F) ^~^ F) \ ^"^ r) /'^^^ F) ^~^ F) \\ 



«--(-^"'4%n^"-(S°'^^S"is)-S4-HE"-£;^E^-4 



1 < i < 15. 



Observe that the operators Aj^ 1 < J < 9, and i?j, 1 < j < p, all commute. 

Lemma 4.4. On the locus % = {(t,s) = 0}, the function f := logr^^(t, s; a,/3; a, 6) satisfies 
the Virasoro constraints 






a^/ 



BjBkf + 



n,- nfc iV 2 . , . , „ , \ 



-^^ = -A,B,f -^-'^ + ^-^{{a,m) + {M ' 
dsY'dtY' p q pq pq^ 

9^/ / " 1\ 2 

9^/ / - 1\ 2 

^4^)^ = -(^. - -)Bkf - -{{a,m) + (6,n) ), 



ty/iere (a, m) = Yl'i=i '^i''^i '^'"■d (/?, n) = Yl^^i f^i^i- 
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Proof. We compute on the locus OC, using ( 4.18[ ) that 



^ + JE E-'" ' 



9Si ^ i>2 \ 1=1 



IS- 



ds 



(0 



1 1=1 "'^i-i/ ^ 1=1 1=1 



1/ ^ ^ 

+ - ( ^ aimi + ^ 6, 



1=1 



1=1 



,{fc) 



+ JE E-i" 



ini 



d 



X 



•^■^l ^ i>2 \ 1=1 



«J0 

(y^i-1 1=1 



+ ^Zt, 



(0 



d 



'-"'i-i 



^ q p 



1=1 



1=1 



d 



ds'; 



(fc) 



+ ^E E-!" 



17^ 



df 

.'-'■^1 ^ i>2 \ 1=1 

1/ " ^ 

+ - ( ^ aimi + ^ bini 



7^ + E*4'>-|r)/ + K^""'"-^'"'''' 



i=l 



+ 



g-1 



.'7-1 



p-1 



■^ ^ \ «=i ' 1=1 ' i=i ' 

^ a;mi + ^ bini 
' 1=1 1=1 

d'^f ruk , 1 



X 



q-l 

1=1 



g-l 



+ ^( - Qi'mj - mq){l - 6j^q) + ^{mi - niq) 



1=1 



p-1 



rUn 



2Y,Mni-np)). 



MJ-I, 



1 = 1 



^p— 1 



Since ^^^^^ (m/ - rrig) = iV - gm^, Yli=i <^l{mi - ^q) = («> "i) and X]f=i Al"-; - np) = (/3, n), 
we obtain 



mj ruk N 2 . x , ,0 \\ 



The proof of the other relations is analogous. 

5. Existence of a PDE for log Pp^q{E; a, b) 



D 



In this section we prove that, under the assumptions ai + . . . , +aq = and bi + - ■ ■ + bp = 0, 
the function log Pp^q{E; a, b), with Pp^q{E;a,b) as defined in (2.2), satisfies a nonlinear PDE, 
the variables being ai, . . . , a^-i, 61, ... , 6p_i and the coordinates of the endpoints of the set 
E, i.e. ci, . . . , C2r- To perform this, we first show that the function / := logT^^(t, s; a, I3] a,b) 

satisfies a system oi ^{p + q){p + q — l) equations on the locus £, containing partial derivatives 
with respect to ai, . . . , aq-i, 61, ... , 6p_i, ai, . . . , ag_i, /3i, . . . , /3p_i and the coordinates of 
the endpoints of the set E. 



X 
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We define the operators 

1=1 
5f = -(l-*..)4 + ^(S-.+g|;). 1<.<P. (5.1) 



1 = 1 



'-' 8 jzi. 8 



1=1 ' z=i ' 

We then have 

•*^=4-?(E"'I: + Ea|:)' isis«. 



^.-'^f-^lE-l; +£«!;)■ '^^^^. 

■^ /=1 ' 1=1 ' 

We also introduce the following notation 

- % + -9/3 = - 1 - 5jp)—- + -22w7f, -db, + -db = - 1 - Sjp)-- + - V 7^, 
^ p dPj p ■f^ d/3i ^ p obj p f^ obi 

o 1 o /. r N 5 1 v^ d ^ 1 _ ,^^,51 v^ d 

9 ^ dak qf-^dai q ^ dak q f^^ dai 

(5.3) 

with 1 < J < p and 1 < k < q. Note the two sets of operators on the first row respectively 
sum to zero as we sum j from 1 to p, and similarly for the operators on the second row, as 
we sum k from 1 to q. With these notations we have the following. 

Theorem 5.1. The function f := logr^ ^{t, s; a, f5; a,b) satisfies the following ^{p + q){p + 
q — 1) equationsn on the locus L 

I -^ p ■' p q pq) Aj I q ■' ■' P q Pq ^ B^ 

= Gf^, l<j<q,l<k<p, 



^{f,g}x=gX{f)-fX{g) 



22 



M. ADLER, P. VAN MOERBEKE, AND D. VANDERSTICHELEN 



{Ai{-d^-d^^)f,AfAif+'^+'^-^} +{Af{-da-d.,)f,AfAif+-^ + 



Gl, l<j<A;<g, (5.4) 



{i?,-(i9,-5,,.)/,i?j^i?,-/+!;^ + ^-^}^^^ + {i?j^(i5,-9,J/,i?,-5,-/+^ + 



!!^_^\ 



p p 



i2Jij£ 






1 < j < k < p, 



where GA, G^j^ and G^^ only depend on f, its derivatives with respect to ai, . . . ,ag-i, 
6i, . . . ,bp-i, and its differentials up to the third order with respect to the operators A'^ , B^ 
and So, evaluated on the locus L. 



Proof. Using equations (4.18) we obtain on the locus % 



d 



ds\ 
d 



(i) 



ds{ 
d 



(j) 



dt\ 
d 



(i) 



log 


fn—ej 


+ek,n 


m+ej 


-ek,n 




^E 




log 


m—ej 


,n-el 


m+ej 


,n+el 




^E 




log 


m,n+ 


5j-efc 


^E 






m,n— 


5j+efc 




E 






'm+el 


,n+ej 



dt\ 



(i) 



T 



m—ei^,n—ej 



X 



X 



X 



X 



A 1 ' rn—e)+el,n 2, - 

^j log -R ^ -('^i - (^k), 

m+ej — ef;,n ^ 

T^ . . . 2 

. , m— e,,n— Efc ^, , , 

^j log -g + "("i + ^k)^ 

m+ej ,n+ei ^ 
T^. . . 9 

-5, log -^-^^ + -(6, -6,), 

rn,n—e'j+el. " 

T^ . . . 9 

-5j log -g ^ + -(afc + &j). 



(5.5) 



4.4 we have on the locus % 

^i log -^^ 



Substituting the first equation in (5.5) into the second equation in (3.7) and using lemma 



A, 



^^Akf + ^{{a,m) + {h,n)) 



T -^ - - ^ 

Similarly, we have 

m,n-ej+efc 

^i i'^ ^:b — ^ — = 



^.-4^/ + 7 + ^ 



TV 



q2 



I + ;| ( (a, m) + (/3, n) 



-(aj-afc). (5.6) 



Bj-\)Bkf^fM^.m)^i,h,n) 



T 



Bi In 



T. 



m,n+ej — ej; 



^.^fc/+^ + 



Tlfc N 



m—eii,n—ej 



B,-h)Akf 



2_ 
pq 



^{{a,m) + {l3,n)) 
{a,m) + (6, n) ) 



P 



{bj - bk), 



T 



r 



m+el,n+ej 
E 



-AkBjf 



mfc 



q pq pq 



Aj In _^ 



m-ej,n-ek 



A, 



vABkf 



^1 



^ '^ {a,m) + (/3,n) 
(a, ?n,) + {b,n) ) 



r 



m+ej ,n+ej. 



-^i^fc/ 



nfe _|_ AT _ _2_| 

q pq pq^ 



{a,m) + {/3,n)) 



(ak + bj), 
{oj + hk). 



(5.7) 
(5.8) 
(5.9) 
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and (5.9)jk- We compute Aj(5.8)kj — Bk(5.9)jk, and we obtain 



Let us denote equations (5.6)-(5.9), with indices chosen as above, by (5.6)jk, (5.7)jk, {5.8)jk 



= Ai 



Bk-^)A,f-l^{{a,m) + {b,n) 



pq 



{aj + hk) 



-A,B,f - ^ - ^ + S - I, ( («, m) + (/3, n) ) p 



Bk 



-AjB.f - ^ 



™j nk _^ N 



q ' pq |((a,m) + (/3,n)) q 



{aj + bk) 



since [^j,-Bfc] = 0. Define a{a,b) = {a,m) + {b,n). As Aj and B^ are first order differential 
operators, we have 



{{Bk-l)A,f + 



2a(a,b) 
pq ' ^J 



.1 >^J ' p ' q pq 






AjBkf 



™i nk j_ N 2a{a,(3) 



g pg 



pq 



IV q I '^■' pq ^ J '^■' p q pq pq J 



B, 2 



A.Bkf 



!ILi _ nfc _|_ W; _ 2a{a,l3) 



pq 



pq 



'' -Aj{aj+bk)+-Bk{aj+bk). 
p q 



(5.10) 



Similarly, we compute, for j ^ k, ^^(5.6)^^ + Aj(5.6]kj and Bk{5.7]jk + Bj(5.7]kj, and we 
obtain 



{{A,-l)Akf + ^^,A,Akf + ^ + ^-§ + '-^}^ 



{{Ak-l)A,f+'^,AkA,f + f + 



mk ,rnj N , 2a{a 



+ 



[A.Akf + ^ + 



an j_ m± _ K _L. 2o-(a,/3) 

q g^ g^ 



A, 4 



(5.11) 



and 



{{B,-l)Bkf+'-^,B,Bkf+-f + f-^ + '-^} 



Bk 



"Hi J. Hk _ N j_ 2a{a,f3) 



BjBkf + ^ + ^ - 4 + 
J ""' p p p^ 



{{Bk-])B,f+'^,BkB,f + f + -f-l, + '^} 
+ ^ r^ 



B, 4 



!ife _|_ !iL _ iV I 2a{a,l3) 



BkBif+'^ + 2 



P 



(5.12) 
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Using (5.2) we compute 

g-l 



p-1 



\ ^1=1 ' 1=1 



dbi 



A 



,9-1 



p-1 



^--(E-a^ + EO 



1=1 



1=1 



AfAif + (D{a,P). 



Similarly we have 



AjBkf = AfB^f + 0{a,P), 



j£ d£ 



B,Bkf = BfB^f + 0{a,f3), 



and 



^.--)^;^/ 
4-^j^fc/ 



1 



^^A^/ + -(9a, --aa)/ + 0(a,/3), 



qy q ' q 



da)f + 0{a,f3), 



pJ 



'3 jBkf 

PJ 



-So - {da, da 

q q 

-■B^-{da^--da)--)Btf+-{d, 

q ' q q) P <J 

-So - (dp^ - -dfs) - -)Aif + -{db, - -db)f + 0(a,/3), 

p ^ p pJ q p 



P 



3 p 



P 



Btf + -{dh,--d,)f + 0{a,P). 



P 



P 



Consequently, on the locus L the equation (5.10) can be written 



{{l^o-{d,,-ld,)-l)A^f + l{d,,-ld,)f+'-^,A^Bif + ^ + r^-^] 



p^Ja^- 



q pq 






Bf 2 2 



A^BH 



^i ^k 



p q pq 



- + - 



p q 



Putting all the terms which do not contain derivatives of / with respect to Qj's or /3j's in the 
left hand side, we obtain 



f^'-{^l' 



1 ^ ^'^ -J J '^ p q pq) Aj' 



{{'-da-da,)B^f,AfB^,f+'^ + '^-^\ 



p q pqi Bf 



where 






2 2 



P 
1^ 1 



■' p q pq/ 



{(^So-^)4/ + 5(afc, -£^6)/ + 



p p 

q q 



1 

— ( 
P 



pq ■' p q pq) A^ 



l\q qJ p ■' q pq ■' p q pq) 
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Similarly, on the locus L, the equations (5.11) and (5.12) can be written 






^i 



G 



B 

jk 



i^l 



a,-af,)Bif.BfBif + ^ + '^-^\ 



p p p^ ) B^ 

L p -^ •' V IP V ' B'^ 



where 






jk 



q\ ■' q q g^ / 

v\q qJ q q Q q q q'^JAf: 

IVo o/-^ a a q'^ -^ q q q^ ) A-^ 



q q 

q qJ q ' q 

-'Bo--)Aff+kda,-^- 

q qJ -^ q q 

p\ ■' " p p p^ J 

p pj p ■' p 

p PJ ■' p p 



Gf,:='-{BfBif+"^^'^-'iy- 



l\v »/ V r> v^ ■' V V v^ ) Bf' 



LVr> %)/ ■' V %> V V %) V } B'r- 



p p p^ 

p p P^ - -J 



U 



In order to obtain a PDE for / = log T^ ^(0; a, b) or for log Pp^q{E, a, b), we need to eliminate 



the partial derivatives of / with respect to ai, . . . , a^-i, /3i, . . . , /3p_i from the equations (5.4) 
in Theorem |5JJ Define 



and Yi = {-di3-d[^,)f\j^, 1 < i < p. 



.1 



X^ = {-da-da,)f\^, l<i<q, 



Note that we have XlLi ^i = ELi ^* = 0, and YlUi ^f = ELi ^i = ^-i- Consequently, 
there are among Af , I < i < q, and Bj" , 1 < j < p, only p + q — 1 linearly independent 
differential operators. Set A^ = S_i - Y11=i ^f and B^ = S_i - ^fj^ Bf . 
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With these notations, the equations (5.4) can be written 



I ^ -^ *=•' p q pqiAf \ k J^ 3 kJ p ^ pq)B^ 

= Gf^, l<j<q,l<k<p 



q q" J Aj- 



G% l<j<k<q, 



i^'^.^X +/R^y..R^R^f + !^ + !^-^l 



{B^,Y„BfBU+'^ + '^--A^,+{Bfn,BfBU+'^ + 



p p^ J B^ ^ -^ ■' p p p 

Gffe, \<3<k<p, 



2JiJ^ 



or 



^Sfc S'fc 



(4) V, - [Bi)^X, - (^44^)4y, + i±^Bicf\BiX, = gf,B, 1 < j < g, 1 < A: < p, 

(yl^^)2Xfc + (^f )'X, - (^A'^cfM^Xk - (^AicfMiX, = g% I < j < k < q, 

(5.13) 
{BffYu + {BifY, - (^-^B'^c%)B'^Yk - [-^Bic%)BiY, = gf„ l<j<k<p 

where 



AB aCtjL f "^i '^fc ^ A aC aC f "^i "^fc -^ 



p q pq ■' ■' Q Q Q 

Uj Uk i\ 
p p p' 



^B ._ r£ r£ f I "-i I "■'^ ^ ('t;i/1^ 



and 



fiAB fiA (~iB 

AB _ jk A _ jk B _ jk /c -, c:\ 

Sjk '■- —AB' 9jk ■- — T' 9jk '■- — B"- ICl.iDJ 

jk jk jk 

We have thus a system of M = 2 (P + 9)(P + 'Z ~ 1) linear equations in the p + q — 2 unknown 
functions Xi, . . . , Xg„i and Yi, . . . , yj,-i and at most ah their first and second order derivatives 
with respect to the independent commuting differential operators Af , 1 < i < q — I, Bj, 
1 < J < p— 1, and S_i. We think at all these quantities as unknowns. At this point, we have 
a system with a smaller number of linear equations then unknowns. The general strategy is 
to keep differentiating the equations and show that at some point we must reach a balance 
between the number of equations and the number of unknowns, leading to the vanishing of a 
determinant at the first point this occurs, which must yield a nontrivial relation. Let Zm be 
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the set of linear equations (5.13), and define 
7K 



Z 



M 



[1 + S_i + At + 



+ A^_i + Sf + 



+ 5p^_i]^Z 



M, 



the set of equations obtained by taking the equations of Zm and all their derivatives up to 



the Kt\i order with respect to the differential operators 23_i, A^, . . . , ^o„i, -Bi 



...,i?^_i. The 



number of equations in Z^ is simply M times the number of monomials of degree K chosen 
from a set oi p + q variables, i.e. 



M 



p + q + K-l 
K 



1 



{p + q) 



{K + p + q-l){K + p + q-2)...{K+l) 
{p + q-2)\ 



The set of equations Z^^ is a set of linear equations in the p + q — 2 unknown functions 
Xi, . . . , Xq^i and Yi, . . . , l^-i and at most all their first, second, . . . , {K + 2)* order deriva- 
tives with respect to the differential operators Af , 1 < i < (7 — 1, -B^, 1 < j < p— 1, and 23„i. 
Let L be the number of unknowns in these equations. Then 

[K +p + q + l){K + p + q) . . .{K + ?,) 



L<{p + q-2) 



p+q+K+2-l 
K + 2 



{p + q-2)- 



{p+q-iy. 



Prom these considerations, it is clear that a sufficient condition to have Card{Z^] 



> L is 



2'" ' "' K\{p 
or, simplifying this expression. 



1, ,(K+p + q-iy. , 

(P + g) \.,./, J >{p + q-2) 



21! 



{K+p + q+l)\ 
{K + 2)\{p + q-iy: 



{x^ - 3x + 4)ii"2 + {-x^ + ^x + i)K - 2x{x'^ - 2x - 1) > 0, 

where we have noted x = p + q. We observe that, with p and q fixed, for K sufficiently large, 
this inequality is satisfied, since x^ — 3x + 4 > 0. Let K* be the smallest value of K such that 
this inequality is satisfied, and note k* the number of equations in Z^f , i.e. 



p + q + K* -I 



k* = -ip + q){p + q-l), j^^ 

and L* the number of unknowns in the set of linear equations Z^^ . Let us note these unknowns 
xi, . . . , XL* ■ Then the system of k* linear equations that we have obtained can be written 

(' \ 

Xi 

[a-ijif]] i<i<k* . =0. 

^<j<L*+i : 

\XL* J 

As k* > L* , we can select the L* + 1 first equations in this system and construct the following 
system 

/I ^ 



|.0«i(/)J l<i<L*+l 



Xi 



0. 



\XL* J 
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But then, necessarily, we have 

det [aij{f)]i<i<L*+i = 0. 

This is a PDE of order {K* + 3) for the function /, with variables ai, . . . , ag_i, bi, . . . , bp-i 
and ci, . . . , C2r- Since / = logr^,.(0; a, b) = logPp^g{E; a, b) + logr? ,^(0; a, b), this yields a 
nonlinear PDE of order K* + 3 in the variables oi, . . . , a^-i, bi, . . . , bp-i and the endpoints of 
E for logPp^q(-E; o, 6), and thus throuhg formula (2.1) for logP^^''"'^'' (a^^ Xi{t) G E), in terms 



of the input logr? ^(0; a, 6), which we think of as known. We thus have proven Theorem 



1.1 



6. Non-intersecting Brownian motions with two starting points and two 

ending points 

In this section we consider a particular situation of the problem studied in the preceding 
sections. Consider A^ non-intersecting Brownian motions xi{t),X2{t), . . . ,XN(t) in M, condi- 
tionned to start at time t = at two different points and to end up at t = 1 in two different 
points, with the coordinates of the starting points and the coordinates of the ending points 
both satisfying a linear condition. In this particular case, all the results of the preceding 
sections hold with p = q = 2. Consequently, by virtue of Theorem |1.1[ we know that the 
probability to find all the particles in a certain set E = IJi=i[c2i-i) C2i] C M at a given time 
< i < 1, satisfies a nonlinear PDE of order 6, the variables being the coordinates of the 
starting and ending points, and the endpoints of the set E. The aim of this section is to 
improve the result of Theorem |1.1| in the particular case when p = q = 2 and to describe this 
PDE more precisely. 

For the sake of clarity, we first recall some notations. Consider A^ non-intersecting Brownian 
motions xi{t), X2{t), . . . , XAr(t) in M, with mi particles leaving from a and m2 particles leaving 
from —a, and rii particles ending in b and n2 particles ending in —b. We denote 

/ (xi(0),...,XAf(0)) = (g, ...,a, -a, . . .,-a ) \ 

ah Xi{t) £ E 



+b,-b 



^^^'-^{allxiit)eE) : 

ni 712 



\ 



m2 

(xi(l),...,xjv(l)) = (&,..., 6, 



/ 



the probability to find all the particles in a set i? C M, at a given time < t < 1. We have 



¥l'Z^{^llx,{t)eE)=P2,2{J--^^E;\l^^^^a, I ^' 




t ' \l l-t 

with the normalized problem defined in (2.2). We deform P2,2{E] a, b) by adding four families 



of extra time variables 

,(1) _ (,m ,(1) ) ,(2) _ (J2) ,(2) ^ 

Jl) _ (M Jl) ^ J2) _ / (2) (2) N 

and the two parameters a, f3 £ C. We have 

P2,2{E;a,b;{t,s),{a,^)) = -^ (6.1) 

'mim2;ni,n2 V''! o, u:, fj, u,u) 
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with T^-^m2;ni,n2(.'tiS;a,(3;a,b) defined in (3.4), and where {t,s) = (t*^^\ t*^^^; s^-'^^ s^^^). 
The function P2^2{E;a,b) is then simply given by P2^2{E;a,b;{t,s),{a,(3)) evaluated 
along the locus L = {{t,s) = (0,0), a = /3 = 0}. We define t he fu nction / 
log r,^^ m2-ni n2 (*' '^' '^' /^' '^' ^)- '^^^ following particular version of Theorem 



5.1 



holds. 



Theorem 6.1. Put X 



2da\L 



and Y 



'2 9/3l£- 



The function f 



log r^j m2-ni n2 (*' '^' '^' /^' '^' ^) Satisfies the following 6 equations on the locus L 



7TT-1 ni 



[AiY,AiBff + 
-[AiY,AiB^f + 

{A^,Y,A^,Bff + f + 
-{A^Y,A^B§f 



mi 722 



ni 
2 

n2 
2 



-I 

4 J Af 
AT. 



{BfX,AfBff+'^ + 
.\B§X,AfBif+'^ + 



ni 
2 



m2 






n2 
2 

2 

n2 
2 



-I = 

A ) Bf 
4 /B,f 



-I 

4: ) Bf 



riAB 
•-^11 1 



riAB 

'-^12 5 



<-^21 ■> 



riAB 
'-^22 5 



{AfX,^f^f/ + ^}^^ - {^fX,vlfyl^/ + ^}^^ = G 



■^12) 



{i?fy,i?fi?2^/ + -}^^^ - {i?fy,i?fi?,^/ + -}^^^ 



where 



A': 



a— /3— 
9a db 



A^ 
^i 



' '^-. + (-IF ^ 



S, 



i? 



i?f 



5 ^5 



9a 



db 



^(--(-)'l 



<-^12' 



1 < J < 2, 
1 < i < 2, 



2V~ ^ ' ^ '' ^a/' "-^ 
and where G^^, G£, and G4p on/y depend on f , its derivatives with respect to ai, . . . ,ag_i, 
6i, . . . ,bp-i, and its differentials up to the third order with respect to the operators Aj, BJ 
and "Bq, evaluated on the locus L. 



The equations in Theorem 6.1 can be written 



i^fX 



A^?X 



A^ log cf2 U^X + Ai log cf2 A^X 



(Bf?Y 



B^ iogcf2 )i?2^y + ( sf iogcf2 )i?f y 



{A 

{BifY 

{AifY-{BffX-[Af\ogct,^)AiY+[B^,\ogct,^)B^,X 

{AifY + {Bffx - (4 logcif)AiY - [Bf log 4^)i?fX 

-{AffY 



9l2, 

- 9ii ) 

- 521 ) 



(6.2) 



iB^2fx + (Af logcti'^AiY + (bI \ogct2^)B^2X = 9^2^ 



-{AifY + {Bifx + [Ai \ogc^i)AiY - (b^ \ogcii)B^X 
where the Cjj's and the gij^s are defined in (5.14) and (5.15) with p = q = 



922 ) 



2. 
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The four differential operators Aj,A2,Bf,B2 are not linearly independent. Indeed, we 
have 

Af+A^ = Bf + B^='B^,, "■ "■ ^ ^r ^r d 



Ai-A{ 



da' 



Bk - Bi 



db' 



Consequently, we will rewrite the six equations (6.2) in terms of the three independent, com- 
muting differential operators ^, ^ and S^i. Before performing this, let us introduce some 
notations. First, we will write 

dF 



da ' dh 



Fb, S_i = Fc, 
I uu 

for a function F. Next, we put 

ri „A r' — r,B r — n^^ , ^ab ^ab , ^ab r< — r.^^ , ^ab ^ab ^ab 

<J-1-— 9l2i ^2 ■— 9l2: (j^3 •— -5ll +521 "512 +522 ' ^^4 •— ^ll +521 ~ 5l2 ~ 522 ' 



1 



AB 



G,:=,{gt{ 



AB 



1 



521 



5^2^ + 52^^), G6:=-^(5fi^+52r+5f2^+52l^) 



and 



Ai:=logCi2, A2:=logc 



B 

125 



„AB„AB 
'-12 ^22 



log 



^AB„AB 
-11 ^^12 
„AB„AB ■ 
^21 "^22 



As := log 



AB„AB 
21 '^12 
„AB„AB ■ 
''11 ^-22 



Ae := log {cff4fcfi'42n 



We then define 



a:=-(-A3c + A66) 



7:= 



A4C + A 



5b) 



/3:=-(A6c-A3b) 
<5:=-(A5, + A4fe), 



and 



d:= -(-A4c + A6a) 
7:=-(A3c + A5a), 



/3:=-(A6c-A4a), 



1 



(Asc + Asa). 



Taking adequate linear combinations of the equations (6.2), and using all these notations, 
we obtain the following equivalent system 

X,e = Gi + 2 Ai.X, + 2 Ai,X„, 
Ytc = G2 + -Asfei; + 2^2cn, 



X,, + Xbb = Ga + 13 X, + aXb + 5Y, + ^Ya, 
Ycc + Yaa = G4 + 8X, + 7X, + ^Y, + AY^, 

^hc — Yac = G5 + —Xc + — X5 — —Yc — —Ya 

„ 7 (5 7 (5 

= Ge + ^X, + -Xb - ^y, - -Ya. 



(6.3) 
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This is a system of linear equations in the variables 

-^a) -^bi -^ci '^act '^hcT ^bb) ^cci ^aj ^b? ^ ci ^acj ^bc: ^aa^ J^cc' 

The coefficients of this linear system depend on logT^_^ ^^.n_^n^{0;a,b) and its partial 
derivatives up to the third order with respect to a, b and the endpoints of E. They 
are highly related. Indeed, there are only twelve non zero independent coefficients 
Q,/3,7, 5, a, /3,7, (5, Ai6, Ale, A2a, A2c- We will now generate new linear equations by apply- 
ing successively the derivatives ^, ^ and 23 _i 
system of 37 equations 



^ to this system. Consider the following 



^1), ^2), ^3), ([6^4), ([6^5), ([6^6), ^6)a, ^6)b, ^G)c, 

(pll)a, (p|l)fe, (pll)e, (p|2)„ (pl2)b, (p|2)e, (|63}3)a, (l63}3)fc, K^3)c, (pl4)„, (p|4)fe, (p|4) 



^5)a, ^5)b, ^5)c, ^6)aa, ^Q)ab, ^6)ac, ^Q)bb, ^^)bc, ^6)c 

(l63|l)5fe - (|63}2),, - (p}5),b, 



^2)aa + (l63|2)ec - ^4.)bc, 

(Pl5)aa + §^5)bb + (P}5)ec + (P|2),b - (p|l)„fe + (pl4)„e " (l63|3)fee, 



(6.3,6) 



abc 



(6.3,6) 



2 

(5 



X (6.3 



aac „ ^ (|6.o 



ifcfe 



l)a6 



2 

7 



X i6.3\l)bc + ^ X (|6.3l2)„a + J X (|6.3 



X (6.3 



)ac + 2 X 



((6.3 



l)a6-(6.3 



(6.3,6) 



bbc 



((6.3 5)b6 + (6.3 



2)a6) 



((6.3,5)bc + (6.3 



2)« 



5) 



7 



((6.3 



l)6c-(6.3 



2)ac)+7^ x(6.3 



2)a6 + | x(|6.3 



2)6c- 



These are linear equations, the variables being all the ffist, second and third order derivatives 
in a, b,cof X and Y, except Xaaa and Yhht- Consequently, there are 36 variables. Constructing 
the vector x := (1, Xa, Xb, Xc, Ya, !&, Yc, . . .)'^ ^ C^^ (the first component being one, followed 
by the 36 variables) , this system of linear equations can be written 



[%•(/)] 



1<«J<37' 



0, 



where aij are differential operators of order less or equal then 6. But then we have necessarily 
that 



det [aijif]] 



1<«J<37 



0. 



This is a PDE for / = logP2,2(-E'; a, b) + logrj^^ mrni n2(0; 0; ci) b). Thus using the structure 
of the 6 equations in Theorem |6.1[ we have obtained a much better result than the one in 
Theorem |1.1[ Indeed, performing in detail the general method described in the proof of 
Theorem |1.1[ one obtains in the case when p = q = 2 a PDE given by a determinant of a 
107 X 107 matrix that is equal to zero. Thus in any particular case, one can do much better 
than the general case in Theorem [TTT] 
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Appendix A. The integral over the full range 

In section 5 we have shown that the function / = logr^^ -(0; a, 6) = log Pp^q{E; a, b) + 
logr? ^(0; a, b) satisfies a nonlinear PDE, and in section 6 we have described this PDE when 
p = q = 2. To obtain a PDE for log Pp^q{E; a, b) it is necessary to evaluate logr? ^(0; a, 6). 
This has been done in the particular case when p = 1 (see the appendix in [5j), but it seems 
harder to evaluate this function when p,q > 2. In this appendix, we conjecture some results 
about the evaluation of logr? ^(0; a, h) when p = q = 2. 



Consider N non-intersecting Brownian motions xi(t),X2(t), • • • ,X]\f{t) in M, with mi par- 
ticles leaving from a and m2 particles leaving from — o, and ni particles ending in b and n2 
particles ending in —b. We suppose rrii = ni and m2 = ^2. We denote 

/ (xi(0),...,XAr(0)) = (a, ...,o,-a, ...,-o) \ 



+b~b 



(all Xi{t) G E) :-- 



ah Xi{t) G E 



V 



nil 



"12 



{xi{l), ... ,xn{1)) = {b, ... ,b,-b,. .. ,-b) 



nil "1,2 / 

the probability to find all the particles in a set ii^ C M, at a given time < t < 1. We have 

P^C,"(an Xiit) €E)= P2,2{E;a,b), 



with the normalized problem defined in (|2.2|), and where 

l2{l-t) 



a :- 



a, 



2t 
l-t 



E 



t{l-t) 



E. 



Notice that ab = 2ab. As shown in (3.3), we have 



P2,2{E-ra,b) 



rlMaM' 



with r^ ^(0; a, b) defined in (3.4). We try to evaluate the function 
/ 



r?.(0;a,6) = det 



V 






0<j<mi 
0<j<mi 

0<i<ni2 
0<j<nii 



f^x^+^e^-^-'^^'^e-'^dx 



0<i<nii 
0<j<ni2 

0<i<ni2 
0<j<ni2 



Define 



We have 



f^i+j[c) 



^i+ie-^+cx^^_ 



d \ «+i 



Aio(c) = V27re^, and fJ-i+jic) = l—j ^o(c) 



Define also the polynomials (Hermite polynomials up to a slight change of variables) 



Pj{x) 



e 2 ^ ^ e 2 . 



dxJ 
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We then have 

((/Xj+j(a + 6))o<i<mi-l \lJ'i+j{0' — b))o<i<mi-l 

\fJ-i+j\—0- + 0))o<i<m2~l \fJ-i+j{ — 0' — b))o<i<m2-l 
0<j<mi-l 0<j<m2-l 

({e^"-''pi+j{d + 6))o<i<mi-i {Pi+jia - b))o<i<m.i-i 

^ 0<j<mi-l __ 0<j<m2-l 

(Pi+j(-a + &))o<j<m2-i [Pi+ji-a - b))o<i<m2~l 

0<j<mi-l 0<j<m2-l 

We will use the following lemma. 

Lemma A.l. Consider the block matrix 

A B 
C D 

with A, D square matrices, and D invertible. Then 

det( ^ ^ j =detDx det {A- BD-^C). 
Proof. Doing row and column perations, we have 
det(^ ^J=det(^ {!) ) ^ ^«t ( D^iC f )=det2?xdet(^-/f^ ^ J 
= det D X det {A - BD-^C) . 



Using this lemma, we have 



r|.(0;S,6) = (27r)^2^e^^("+^) e-^^i'^^detD x det (^ - ^D-^C), 
where 



-1' 

0<j<mi-l " ' 0<j<m2-l 



D 



0<jr<mi-l 0<j<m2-l 

and it is well-known that det D = 01=0^ ^ - ^^^ ^^ note 

^ ^ ^4.^ _J^'{^^ j^ ± f e^dz ^ l_ r e^dz 

U ^- 2vrz7r^^^,^-2(,_4a6) 27rz Jr,^^^^ z™^ (z - 8a6) ' 

(A.2) 

where r„ ^-^ denotes a closed contour containing and 4o6 in the complex plane. Computer 
observations point out that for A, B, C, D as in (|A.l) we have 



det(^-5D-iC) =detfpi+,(a + 6)A + Pij(a,6)jo<,<„^^_i, (A.3) 

0<j<mi-l 
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where Pij{d,b) is a polynomial in a,b such that Pij{d,b) = Pj,i{h, a), P ofi{d, b) = 0, and the 



order of Pjj(a, b) is 2(mi — l)+i+j when i+j > 0. We will develop (A. 3) in the large m2-liniit 




Figure 1. Non-intersecting Brownian motions in the large ?Ti-2-limit, with mi fixed 



Let us consider the large m2-limit, keeping mi fixed, see Figure [TJ If mi = 0, for large m,2 
the mean density of brownian particles at any time < t < 1 is supported by an interval with 
endpoints given by iby^2m2t(l — t) — a(l — t) — bt. When m,i is fixed but not necessarily zero, 
the non-intersecting nature of the cloud of rn-2 particles implies that the largest one will again 
reach a height of about \/^ — ^2~ Sit t = g- We will consider the following scaling given in 
[2] for the starting and the target points 

B 




'-2 "^2 

With this scaling, the m^i wanderers will interact with the bulk of m,2 particles (m2 very 
large) , upon considering regions close to x = y/^ — ^^ and t = ^, namely at space-time 
positions (x, t) which scale like 



t 



1 1 

2^2 



m 



1/3' 



X 




m2 
2 



+ 



e 



+ 



.2 . 2^/2m^/'^ 4^2 

We suppose A < B. Under this scaling, the quantity 

A-B 



m 



1/6 



iB-A) + 



1 {A + B)t 
4^2 



m, 



1/6 



Sab = m,2 1 + 



m,. 



1/3 



AB 

2/3 



m, 



•'2 '"-2 

is strictly less than ?n,2, for ?Ti2 large enough. Consequently, by Cauchy's theorem, the contour 



ro,8a6 in (A.2) can be taken to be a circle centered at the origin and of radius m,2- We will 
follow [2j to obtain an assymptotic expansion for X. Making the change of variable z = um.2 
in the integral definig X we have 



X = {Sab^'^ 



m> 



— m2 



2TTi 



r,m2F(u) 



\u\=i -a _ 1 _ (^ _ B)m2 ^'^ + ABm. 



-2/3 



du. 



(A.5) 
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where 



with 



F{u) := n - Inn = 1 + -{u - if + 0{{u - if), 



R{F{u)) =3fi(u)-ln(|u|). 



(A.6) 



The stationary points of the function F{u) are solution of the equation F'{u) = 0, and thus 
there is one stationary point at n = 1. We can deform the path |t(| = 1 into 75 = {l + iy\ — 5 < 
y < ^} plus a circle segment 7' centered at the origin and joining the extremities of 75. It 
follows from (A.6) that 75 is tangent to the steep descent path through u = 1. We choose 



-2/5 



6 = 777-2 • Then the contribution to the integral coming from 75 is given by 

m2F{u) _my^e^'^ fl+ira^ 

2/5 



15U-1-{A- B)m2 ^^^ + ABm^ ^^^ 
Making the change of variable w = — ii 

r Qm2F{u) 

J-ys u-l-{A- B)m~^^^ + ABm^^^^ 



du 



-rrir, 



(A-B) 



l—im^ 



e^(«-l)'d77(l + 0(777-^/^)). 



irn^ 



{u — 1), we obtain 



du 



iV2' 



m 



-1/6, 



,rra2 



As 



we have 



+00 
1 1/10 



I 



^m2F{u) 



■ysu-l-{A- B)m^ ^^^ + ABm^ ^^^ 



{A-B 
dio = o (7772 ) , 

du - 



1 1/10 
75 ™2 



1/5n 



e-" da; (1 + 0(777^'/')) 

1 1/10 ^ \ z // 



75 '"a 



7V27r777. 



^1/6, 



,m2 



;i + 0(7772 ^/^)) 



l/5^ 



U 



(A-B) 

Let us now evaluate the contribution to the integral coming from 7'. Along 7', we have 
= Vl + b'^e'^ with I cos 01 < 



2/5 



and (5 = 7770 , and thus 



|gm2F(«)| ^ gm25R(F(«)) ^ ^022 (v^I+P cose-| ln(l+52)) ^ ^^2^-^ ln(l+<52) 

It follows that 



Along 7', we also have 



g -m2|gm2F(M)| _ Qi^-\^<^2^ 



< 



7i - 1 - (^ - 5)7772 + ABm^ 

and thus 

1 



-1/3 



l^/TTF - 1 - 1^ - fi|7772 "'"^ - l^fi|777; 



-2/3| 



77 - 1 - (^ - S)7772 ^ + ^^7772 



< 



-777- 



1/3 



A-B 



[1 + 0{m-'^')) 



l/3s 
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1/5, 



It follows that the contribution to the integral in (A. 5) from 7' is of order 0{e 

then the main contribution coming from 75, for some < c < 2- Consequently we have 



smaller 



e^ dz 



1 



"^'"^'^'^'i^riA-B) 



l + 0(m-^/^)). 



l/5x 



It follows that in the large ?7i2-limit 

X = J^'^_ r exp (ms + {A- B)mf^ - -(yl^ + B^)m)l^ + -AB{A - B)) (l + Oim^^^^)) . 



/2tt{A-B) """^" \ '^ ' "'' ' ^ 2' ' ' 2 

Consequently, in the large ?7i2-limit, for mi fixed, T?^(a, 6) is expected to behave like 



??i]^+Tn2 m\-\-ra'2 , 



mi —1 



Tl^id,b)^i2n)-^e-^^'^'^\-'"^^'\ n ^0( n i!)^™^- 

j=0 



m2-l 

n 

j=0 
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